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Abstract. We demonstrate generation and frequency doubling of unit charge vortices 
in a linear astigmatic resonator. Topological instability of the double charge harmonic 
vortices leads to well separated vortex cores that are shown to rotate, and become 
anisotropic, as the resonator is tuned across resonance. 
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Optical vortices are topological objects whose transformation properties under 
propagation in linear and nonlinear optical media have been the subject of much recent 
work[lJ. The vortex charge of a beam, defined as the closed loop contour integral of 
the wave phase modulo 27r, is generally a conserved quantity under linear propagation 
in free space. Optical vortices occur naturally in speckle fields^, and can be generated 
in a controlled fashion using diffraction from holographic plates [3J Ej. They can also 
be generated in lasers jSJ E] and cavities with nonlinear elements[7], while high order 
vortex modes have been observed in active cavities with field rotating elements [HI E|- 
Astigmatic optical elements as well as nonlinear wave interactions can be used to change 
the vortex charge of a beam. For example in the weak pump depletion regime of 
second harmonic generation the amplitude of the envelope of the harmonic field at 
frequency co>2 can be written as A2 ~ A\. Thus an input field with charge m of the form 
A\ ~ e m ^ generates an output field A2 ~ e l2m 4> w jth twice the charge. This effect has 
been demonstrated experimentally by several groups using optical vortices created by 
diffraction from a hologram, and then allowing the beam to pass through a frequency 
doubling crystal|HIl HQ 112 13 d ■ 

In this work we describe a different approach to the generation of vortices in second 
harmonic generation that is based on frequency doubling of resonator modes with a 
vortical structure. Consider an empty resonator with the pump beam mode matched to 
the lowest order transverse resonator mode (TEMoo mode), which has a slowly varying 
amplitude at the cavity waist given by u 00 ~ exp(— r 2 /w 2 ), with f = xx + yy and w c 
the cavity waist. By changing the cavity tuning, and slightly tilting and displacing 
the pump beam, we can couple to higher order transverse modes of the cavity. By 
appropriate alignment of the pump beam it is possible to couple to a single higher order 
transverse mode, such as Ui ~ xexp(— r 2 /w 2 ) or u 01 ~ yexp(— r 2 /w 2 ) which have edge 
dislocations. Following the cavity by an astigmatic mode-converter [T5] the Hermite- 
Gauss modes can be efficiently converted into azimuthally symmetric Laguerre- Gauss 
modes with non-zero vortex charge, as was demonstrated by Snadden et al.|16j. 

As we show here it is also possible to generate a vortex mode directly, without 
using an astigmatic mode converter, by aligning the pump beam to give the desired 
superposition of wio and woi modes. Let the pump beam be a displaced and tilted 
Gaussian. Ignoring unimportant constant amplitudes as well as any overall phase 
we have u p = exp (— \f — r p \ 2 /w 2 ) exp [iq ■ r). Here w p is the pump beam waist, r p = 
x p x + y p y is the transverse displacement of the pump beam, and q = q x x + q y y is the 
transverse wavevector that is proportional to the pump beam tilt in the x, y plane. The 
lowest order odd cavity mode that the pump couples to can be written as 

u = o w a(v - u w )uio + o 01 a(v - v 01 )u 01 (1) 

where o mn ~ / dxdy u*u mn is an overlap integral, and a(u — u mn ) is a complex 
coefficient that depends on the difference between the pump frequency v and the 
resonant frequency of the mode v mn . The amplitude and phase of the overlap coefficients 
can be independently varied by adjusting the pump beam^T]. Performing the integrals 
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Figure 1. Experimental setup. 



we find oin = h(2x p /w p — ik x w p ) and oni = h(2y p / Wp — ikyWp) , where ft, is an unimportant 
common factor. By adjusting the pump beam we can obtain oqi = Oioexp (z7r/2) which 
results in vortex generation when a(V — u w ) = a(v — i/ i). 

Vortices were generated in this way using the experimental setup shown in Fig. ^ 
A Ti: Sapphire laser at 858 nm generates a continuous wave fundamental beam with a 
power of up to 300 mW incident on the freqency doubling cavity. The beam is mode 
matched to a linear cavity with two R = 25 mm end mirrors that contains a 1 cm long 
a-cut KNbC>3 crystal with anti-reflection coated ends. The input and output mirrors 
had T 858nm = 4.8%, 0.04% and T 42 9 nm = 92% so that only the fundamental field was 
resonant in the cavity. Phase matching was controlled by varying the temperature of the 
crystal. With the distance between the mirrors set to L cf = R+ ^L c = 30.7 mm for 
confocal operation (n c , L c are the crystal refractive index and length), and the crystal 
temperature tuned for large phase mismatch so no harmonic beam was generated, a 
cavity finesse of about 80 was measured. This agrees well with the theoretical value 
of F = 84 that was calculated using measured values of the crystal losses. With the 
crystal temperature tuned for optimum phase matching up to about 60 mW of 429 nm 
light was generated in a TEM o mode. 

The cavity length was then reduced by 1.63 mm which resulted in the appearance 
of higher order transverse modes in the cavity transmission spectrum as seen in Fig. 
El Using the measured free spectral range as a scaling parameter the theoretically 
calculated frequencies of the first few higher order (q, m, n) modes have been indicated 
in the figure. The observed resonance frequencies agree to within a few percent with the 
calculated values. The cavity was locked to a (q, 1, 0) resonance using the rf sideband 
technique PH] which resulted in stable generation of a unit charge vortex mode as seen 
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Figure 2. Scan of resonator transmission with theoretical positions of the modes 
marked by arrows, q denotes the axial mode index. The upper line is the linear ramp 
voltage applied to the cavity piezo,and the inset shows the far-field structure of the 
fundamental field at the odd mode resonance. 



in the inset of Fig. |21 

The phase structure of the harmonic field was observed by interference with a 
TEMoo beam generated in a second doubling cavity. As expected the second harmonic 
beam contained a doubly charged vortex as seen in Fig. El The detailed structure of 
the beam had a sensitive dependance on resonator alignment. Careful alignment of the 
crystal position resulted in observation of a doubly charged core region, although there 
was an apparent tendency for the vortices to repel each other so that a small vortex 
separation remained as seen in Fig. El We attribute the splitting to the topological 
instability of m > 1 vortices [TH]. Adjustments to the crystal and/or pump beam 
alignment resulted in the core splitting into two well separated singly charged vortices, 
with a controllable relative orientation and separation. It was also found that when 
the resonator was aligned so that the cores were well separated, the relative orientation 
angle of the cores rotated in a repeatable fashion as the resonator was tuned across the 
(q, 1, 0) resonance, as seen in Fig. HI 

The observations of vortex rotation can be explained by taking account of the 
cavity astigmatism due to crystal birefringence. Let the fundamental beam propagate 
along z (a-axis of KNbOs) and be polarized along y (b-axis). Beams propagating in 
the x — z plane correspond to ordinary polarized rays with an index n b = 2.279 at 
858 nm and the temperature set for phase matching. We will label the x — z plane 
modes with the first transverse index m (the x axis is horizontal in the figures). Beams 
propagating in the y — z plane correspond to extraordinary polarized rays with an index 
n(6) = n b (l + tan 2 ^) 1 / 2 /(l + (n b /n a ) 2 tan 2 #) 1//2 where 6 is the angle of the ray with 
respect to the y axis and n a = 2.238 We will label the y — z plane modes with the 
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Figure 3. Near field images of the harmonic charge 2 vortex. The central frame is 
overexposed to reveal the vortex splitting, as verified by the far-field interferogram on 
the right. 



second transverse index n. The effective crystal thickness for x — z plane modes is L c /nt, 
while the effective thickness for y — z plane modes is L c /n, where n = n 2 a /nb is the 
effective index |2*U]. Since n < n b the crystal is effectively longer in the y — z plane 
and these modes have lower resonance frequencies. A short calculation shows that the 
frequency splitting can be written as 




(2) 



^qmO ^qOn 

I 

n b , 

where L = R + [{n 2 c — l)/n c ]L c + 5 with S = L — L c f the change in cavity length 
from confocality. For our exprimental parameters the shift given by Eq. (0) is about 
15 MHz, which is a non-negligible fraction of the cavity resonance which has a FWHM 
of 27 = 27r x 42 MHz. Indeed, close inspection of Fig. |21 reveals that the odd order 
transverse modes are somewhat wider than the lowest order modes, in agreement with 
Eq. (j2J). The modal coefficient a that appears in Eq. can be written as 

a(u - u mn ) = [l + (4F 2 /vr 2 ) 2 sin 2 (7r(^ _ u mn ) / {2 1 F))]~ 1 ' 2 e^'^-\ (3) 

where Xmn = fsm[iv(u - u mn )/(Fj))/{l - f cos[ir(v - u mn )/(Fj))}, and f = 1-n/F. 
For our experimental parameters the relative phase shift when the resonator is tuned to 
the midpoint between the (10) and (01) resonances is tan" 1 xio — tan -1 X01 ~ 40°. 

The simulations of vortex rotation shown in Fig. HJwere obtained by calculating u 
from Eq. (0), squaring the field, and adding a small constant offset which resulted in 
splitting of the vortices. The phase difference between 010 and 001 was set to 7r/2 — 40° to 
give the correct 7r/2 phase shift for v at the midpoint between the modes, after which the 
frequency v was tuned across the resonance. We observe qualitative agreement between 
the observed rotation, and the calculated results. The reason for the calculated rotation 
being somewhat smaller than that observed may be that the calculations correspond 
to the near field while the observations shown in Fig. 0]were recorded in the far field. 
Additional observations revealed that the rotation angle was noticeably smaller in the 
near than in the far field. This is fully consistent with the propagation induced rotation 
of a vortex pairJUIj. It should be emphasized however that the observed dependence 
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Figure 4. Rotation of harmonic vortices observed in far-field (top row) and calculated 
from Eq. Q (bottom row). The columns are labeled with the pump beam frequency 
which increases from left to right. 



of the rotation angle on the cavity tuning is not a propagation effect, but rather due 
to the interference of nondegenerate cavity modes with amplitudes given by Eq. ©• 
We also note that the shape of the vortex cores becomes elliptical on either side of the 
resonance. This may be due to the unequal changes in amplitude and phase of a which 
implies that we are generating and doubling anisotropic vortices |2T|. 

In conclusion we have observed creation, doubling, and splitting of vortices in a 
second harmonic generating resonator. Tuning of the resonator across a pair of non- 
degenerate transverse modes results in rotation of the harmonic vortex pattern. 
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